Magnetic exchange interaction induced by a Josephson current 
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We show that a Josephson current flowing through a ferromagnet-normal-metal-ferromagnet tri- 
layer connected to two superconducting electrodes induces an equilibrium exchange interaction be- 
tween the magnetic moments of the ferromagnetic layers. The sign and magnitude of the interaction 
can be controlled by the phase difference between the order parameters of the two superconductors. 
We present a general framework to calculate the Josephson current induced magnetic exchange in- 
teraction in terms of the scattering matrices of the different layers. The effect should be observable 
as the periodic switching of the relative orientation of the magnetic moments of the ferromagnetic 
layers in the ac Josephson effect. 

PACS numbers: 75.70.Pa., 75.30.ds, 73.40.-c, 75.70.-i 






a 



X3 

o 



> 

oo 
in 

o 



cd 



I 

o 
o 



X 



I. INTRODUCTION 

Despite their apparent simplicity, ferromagnet- 
normal-metal-ferromagnet trilayers exhibit many inter- 
esting properties. One example is the equilibrium ex- 
change interaction between the magnetic moments of the 
two ferromagnets, which is mediated by the coherent elec- 
tron motion in the normal metal spacer layer .H Depend- 
ing on the thickness of the spacer layer, this interaction 
may favor parallel or antiparallel alignment of the mag- 
netic moin£p.ts,Ba or, in some cases, even perpendicular 
alignment .Eru In addition, a torque may be exerted on the 
magnetic moments when an electrical current is passed 
through the trilayer. For this nonequilibrium magnetic 
torque, the preferred magnetic configuration (parallel or 
antiparaJld) was found to depend on the sign of the 
current ,Bij so that a reversal of the current switches the 
magnetic moment of the ferromagnets from parallel to 
anti-parallel. This reversal can be observed by a mea- 
surement of the conductance, which is larger in the par- 
allel configuration than in the antiparallel one (this is 
known as giant magneto resistancetl) . 

A unified description of equilibrium and nonequilib- 
rium torques can be given using the concept of spin 
current. In the case of a ferromagnet-normal-metal- 
ferromagnet (FNF) trilayer, this description was intro- 
duced by Slonczewski as an alternatiiie way to calculate 
the equilibrium exchange interactionJiJ (The "standard" 
way to calculate the exchange interaction involves a com- 
putation of the derivative of the free energy to the angle 
between the two magnetic moments.) When electrons 
scatter from a spin-dependent potential, as is appropri- 
ate for a mean-field description of ferromagnetism, the 
spin current carried by the conduction electrons need not 
be conserved. Since the total spin of the system (i.e., the 
combined moment of the conduction electrons and the 
other electrons responsible for the magnetic moment in 
the ferromagnetic layers) is conserved, the lost spin cur- 
rent must have been transferred to the magnetic moment 
of the ferromagnet, which means that a torque is exerted 
on the moments of the ferromagnets. In this way, the 



equilibrium exchange interaction is seen to follow from an 
equilibrium, spin current flowing from one ferromagnet to 
the otherjlj much like the equilibrium (jiersistent) cur- 
rent that exist in mesoscopic metal ringsjl^l whereas the 
nonequilibrium torque arises from the non-conservation 
of spin ciirrcnts flowing in conjunction with the electrical 
current £1 

The nonequilibrium torque is changed when the FNF 
trilayer is coupled to one superconductor (S) contact and 
one normal-metal (N) contact, instead of to two normal- 
metal contacts. The main difference between N and S 
contacts is that the former can carry both spin and charge 
currents, while the latter can only carry a charge current 
for voltages below the superconducting gap Aq. In a pre- 
vious publication, we showed that this restriction gives 
rise to a nonequilibrium torque that, depending on the 
direction of the electrical current, can lead to the switch- 
ing of the magnetic moments to a perpendicular co, 
figuration, rather than a parallel or antiparallel one 
The equilibrium torque (i.e., the magnetic exchange in- 
teraction), however, is not qualitatively affected by the 
presence of the one superconducting contact .£3 

In this paper, we consider an FNF junction with two 
superconducting contacts. For this system, a macro- 
scopic supercurrent may flow through the junction al- 
ready in equilibrium, the magnitude and sign of the 
current depending on the phases of the order parame- 
ters of the two superconductors. As ferromagnets break 
time-reversal symmetry, they are expected to suppress 
the Josephson effect. However for sufficiently thin or 
weak ferromagnetic layers (for instance a Cui_a:Nia; al- 
loytJ with X > 0.44) , the Josephson effect may sur- 
vive. Magnetic Josephson junctions with one ferromag- 
netic layer have received considerable attentipn because 
of the possibility of 7r-junction behavior Jl3^ 
been observed experimentally only recently, 
son junctions with two magnetic layers were studied in 
Refs. po| , pi| , where it was shown that the supercurrent for 
antiparallel alignment of the magnetic moments can be 
larger than for parallel alignment. 

Here, we consider the exchange interaction in an FNF 
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junction with two superconducting contacts. We find 
that the equihbrium exchange interaction is deeply af- 
fected by the presence of the superconductors. By 
the same mechanism by which the supercurrent de- 
pends on. the, relative orientation of the two magnetic 
momentSjoO the exchange interaction depends on the 
phase difference between the two superconducting order 
parameters. As a result, the supercurrent controls the ex- 
change interaction between the two magnetic moments. 
In contrast to the usual magnetic exchange interaction, 
which involves contributions from the whole energy band, 
the Josephson current induced magnetic exchange inter- 
action is carried only by states with an energy within a 
distance of order Aq from the Fermi energy. At a first 
glance, one is tempted to consider the Josephson current 
induced torque as the direct analogue of the nonequilib- 
rium current-induced torque that exists for normal metal 
contacts. However, as we'll show in the remainder of this 
paper, that is not the case: Apart from its magnitude, 
the Josephson current induced torque has most of the 
features of the standard equilibrium magnetic exchange 
interaction. 

This article is devoted to the study of the Josephson 
current induced magnetic exchange interaction and to its 
consequences on the dynamics of the magnetic moments. 
It is organized as follows: In Section II, we present the 
concept of spin current and discuss the differences be- 
tween equilibrium and nonequilibrium torque. We then 
focus on the case of a Josephson junction (both electrodes 
superconducting) . Section III contains a brief review of 
the scattering approach, that allows for practical calcu- 
lation of the torques discussed in Section II. We are then 
ready to discuss, in Section IV, the magnetic exchange in- 
teraction in the Josephson junction, using various models 
for the scattering matrices of the normal and ferromag- 
netic layers involved. Finally, the effect of the torque 
on the dynamics of the magnetic moments is briefly dis- 
cussed in Section V. 
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FIG. 1. Schematic of the trilayer system considered in this 
article. Fa and Fb are the ferromagnetic layers and rfia and mj, 
their respective magnetic moments. The ferromagnetic layers 
are separated by a normal spacer A'' and are connected to two 
reservoirs that can be either normal [N) or superconducting 
(5). The numbers 1, 2, 3, 4 stand for ideal fictitious leads 
that have been added for technical convenience. 



II. SPIN CURRENT AND SPIN TORQUE 

The system we consider is shown in Fig. y. It 
consists of a ferromagnet-normal-metal-ferromagnet tri- 
layer, connected to (possibly superconducting) electrodes 
on the right and the left. The two ferromagnetic layers 
are labeled Fa and Fb, the normal metal spacer layer is 
labeled N . The magnetic moments of Fa and Fi, point 
in the direction of unit vectors rha and ?fif,, respectively. 
The angle between rha and mb is 9. We assume that rhh 
points in the z-direction. For technical convenience, we 
have added pieces of ideal lead (labeled 1, 2, 3, and 4) 
between the layers F^, N , and F^, and the reservoirs. 

The conduptjon electrons are described by an effective 
Hamiltonian,E3 



H,s^ fdr J2 *aW^"/3*/3(r) 



(1) 



a,/3=T,i 

df [A(f)* j (f)*| (f) + A* (f)*x (^)*T (^) 



where '^a{r) creates an electron with spin a and A(r) is 
the superconducting gap. In the normal regions, A{f) — 
0. Finally, the 2x2 matrix 

H = -{h^/2m)\/'^ + V{r) - Ep 

contains kinetic, potential, and Fermi energy. The po- 
tential V{r) represents the spin- independent scattering 
from impurities, as well as the spin-dependent effect of 
the local exchange field inside the ferromagnets. For sim- 
plicity, we assume that the local exchange field is always 
parallel to the total magnetization of the layer (which 
corresponds to the neglect of spin- flip scattering). Thus 
we have 



y(f) 



Vr 



maj (r) 

V;„in(r) 



(2) 



Here, maj (min) stands for majority (minority) and 
exp[— «crj,6'(r)/2] is a rotation matrix rotating from the 
external reference frame to the direction of the local ex- 



change field. Outside the ferromagnets, V^ 



K, 



For the system under consideration, the angle f?(r) = 9 
inside Fa and zero elsewhere. 

An expression for the spin current is obtained by writ- 
ing down the conservation equation for the spin density 
r]{f), 

^>1 = |E*"(^1^"'3*/5(r1, (3) 

a/3 

<? = {o'x, (Ty, cTz)'^ being the vector of pauli matrices. The 
time evolution of spin density reads, 

= -V- 7 [r] + '-{^Hr)[V{f),d]^{r)), (4) 



with the spin current density tensor j defined as, 

j (r) = -^(*^(f)VCT*(f)- V*Vf)(7*(r)). (5) 

The spin current density has one index in spin space 
and one in real space while the brackets (. . .) stands for 
the quantum mechanical expectation value. Equation (0) 
shows that, unlike charge current, spin current is not con- 
served inside the ferromagnets. In fact, the current in- 
duced torque and the magnetic exchange interaction fol- 
low from the spin current lost by the conduction electrons 
inside the ferromagnets. The current induced torque fol- 
lows from the non-conservaton of non-equilibrium spin 
current, while the magnetic exchange interaction follows 
from the non-conservation of the equilibrium spin current 
between Fa and Fb- The total torque Ta, fb (i.e., the sum 
of equilibrium and non-equilibrium contributions) on the 
layers Fa and Fb is found as 



fa — Ji — J2 , fb — J3 — J4 



(6) 



where Ji is the total spin current that flows in the x 
direction in region i (i = I, 2, 3, 4), 



A 



dydz j^{x,y,z) 



X <E i. 



(7) 



We now focus on the torque fa on the magnetic mo- 
ment of layer Fa. In addition to the unit vector rha that 
points along the magnetization direction of Fa , we intro- 
duce the unit vectors u = rha x "T-ft/l'^ia x 'T^bl, which 
points normal to the plane spanned by rha and to;,, and 
w = rha X v, which lies in the plane spanned by rha and 
rhb, but points perpendicular to rha, see Fig. 0. For the 
configuration of Fig. y, {) is the unit vector in the y di- 
rection and the plane spanned by rha and rhb is the xz 
plane. From the observation that 



dV 



(8) 



combined with Eq. (0) , we find that the out-of-plane com- 
ponent of the torque r^ = Tq • u is equal to the derivative 
of the energy E = {H) of the trilayer to the angle 9. 
No such simple result can be found for the torque in the 
w-direction (the component of the torque in the plane 
spanned by TOq and rhb). Hence, we find for the total 
torque acting on the magnetic moment of Fa 






(9a) 
(9b) 



T„ = —-, 



dr (*t(^(y^^^. _ ^^^.^(^ . ^))^(f)). (9c) 



fS-Fa 



Since the spin current in the direction of TOq is always 
conserved in the absence of spin-flip scattering, there is 
no component of the torque along rha, cf. Eq. (0). Al- 
though the above derivation may seem a little specific. 




FIG. 2. Schematic of the out of plane component of 
the torque t^v (left side) and the in-plane component of the 
torque t^w (right side). The equilibrium exchange interac- 
tion only has an out of plane component while the non equi- 
librium torque is mainly in plane. 

Eq. ^ can be shown to hold as well in the presence of 
two-body interactions. 

As seen from Eq. (g), the total torque consists of two 
contributions, which are sketched in Fig. 0. The in-plane 
torque r™ directly pushes rria towards or away from rhb- 
It is the main component of the nonequilibrium torque, 
which is discussed in Refs. BbI In equilibrium, however, 
no spin current can flow outside the trilayer, so that, by 
Eq. d). 



=:cqu 7 7 ^QU 



(10) 



Combined with the requirement that fa (fb) is perpen- 
dicular to rha {iTib), this relation impli^ that, in equilib- 
rium, the in-plane torque r™ vanishes. t3 The out-of-plane 
torque r^ causes a precession of one magnetic moment 
around the other one. This is similar to the Larmor pre- 
cession of the moments in a (possibly 0-dependent) mag- 
netic field. In the presence of dissipation, the system will 
then relax to the lowest energy configuration, where the 
energy is minimal and, hence, by Eq. (9b), the torque 
zero. 

According to Eqs. (|l^) and (^), there are two, equiv- 
alent, ways to calculate the equilibrium torque f^'^^: As 
the equilibrium spin current J2 flowing from Fa to Fb, 
or as the derivative of the ground state energy E to the 
angle 9 between rha and rhb- (Note that there is a direct 
analogy between the equilibrium spin current flowing in- 
side the trilayer with an angle 9 between the two mag- 
netic moments and the persistent current in a mesoscopic 
ring in presence of an Aharonov-Bohm flux.) In the re- 
mainder of this paper, we concentrate on the equilibrium 
torque in the case where both left and right electrodes 
are superconducting. 

The effective Hamiltonian Heff can be diagonalized by 
the Bogoliubov transformation. 



*a(r) = ^ Ua.eif)-i^ + <,e(^)7l, 



(11) 



e>0 



where the fermion operator 7J creates an excitation at 
energy e and u {v) is the electron (hole) component of 
the solutjiep to the Bogoliubov-De Gennes equation at 
energy ejHS 



H 

-iA*cr,, 



lAcTj, 

-H* 



(12) 



In principle, the superconducting order parameter A(r) 
has to be calculated self-consistently, 



A(r) = -^|5(r)|^[l-2/(e)]«V., 



yU, 



(13) 



e>0 



where g{r) is the BCS interaction constant, which is finite 
inside the superconductor and drops abruptly to zero out- 
side S, and /(e) is the Fermi function. Here, we adopt the 
simple model that A(r) has its bulk value A = Age*"^/^ 
(A = Aoe^*'^/'^) inside the left (right) superconducting 
reservoirs, while A = in the normal layers. At the 
normal-metal-superconductor interface, A(r) can be ap- 
proximated by a step function. This approximation, dis- 
cussed in Ref. ESl is valid for the quasi-one dimensional 
geometry we consider here. Up to terms that are inde- 
pendent of 9 and 0, iJoff reads 



^«ff = -?E^+E^ 



ihe- 



(14) 



e>0 e>0 



Combining Eqs. (^ and ([Tl|), one obtains for the spin 
current 



Im Y^ {v)Vdu - v'^VavA {^l^i, - i) . (15) 



e>0 



2m 



Equation (113) was used to calculate the non- 
equilibrium torque for an FNF trilayer with one super- 
conducting contact in Ref. |l^. In that case, only eigen- 
states of the Bogoliubov-de Gennes equation with en- 
ergy close to the Fermi level were involved, so that an 
analog of the Landauer formula could be derived to find 
the noncquilibrium current-induced torque. The equi- 
librium torque, however, is of a very different nature 
since it involves all states in the entire conduction band. 
This limits theoretical approaches to analytical calcu- 
lations for simjJe model systemil or ab-initio numeri- 
cal simulations.B This difficulty does not occur for the 
Josephson current induced torque. 

In a Josephson junction, the equilibrium current / at 
finite temperature is given by the derivative of the free 
energy F of the junction to th£-|phase difference </> be- 
tween the two superconductors ,t3 



/ 



2edF 



(16) 



This equation is very similar to the equation for the equi- 
librium spin current between the two ferromagnetic lay- 
ers, J2 = 



equ 



where Eq. (||) gives 



dF 



Combining these last two equations one finds that, 



(17) 



dl _ 2e ^T^i" 
'do ~ ¥ d(h ■ 



(18) 



In other words, a ^-dependence of the supercurrent im- 
plies a (/)-dependence of the equilibrium torque. This is a 
very suggestive result, as it was recently predicted that 
the Josephson current should be very sensitiue-to the an- 
gle between the two magnetic momentsHj^Ell In what 
follows, we refrain from calculating the full equilibrium 
torque but concentrate on it 0-dependent part, the ratio- 
nale being that it is precisely the (/)-dependent part that 
can be viewed as the supercurrent induced torque. As 
we shall see, the (/)-dependent part of the torque only has 
contributions from energies within Aq of the Fermi en- 
ergy and can, therefore, be calculated from the scattering 
properties of the junction at and near the Fermi level. 



III. SCATTERING MATRIX FORMALISM 

A general review of the scattering matrix formalism 
can be found in Ref. ^, while the particular application 
to the calculation of spin currents in FNF trilayers is dis- 
cussed in Ref. H. Below we give a brief compilation of 
the necessary formulas for the calculation of the equilib- 
rium torque in an FNF trilayer with two superconducting 
contacts. 

The trilayer is bounded in the y and z directions, so 
that the corresponding degrees of freedom are quantized 
and give rise to A'^ch propagating modes at the Fermi 
level, with A^ch ^ A/X'^p, A being the cross section of 
the junction and Xp the Fermi wave length. Each trans- 
verse mode appears as a left moving mode and as a right 
moving mode, and with components for particle/hole 
and spin degrees of freedom. We expand the solution of 
the BdG equation in terms of these modes and describe 
wavefunctions in terms of the 4A'ch-component vectors 
^j which are the projection of the wave functions 
(u, v) on the left (right) going modes in the ideal lead 
i (i = 1, 2, 3, 4, see Fig. 0). The layers F^, ft, and A^ are 
characterized by 8 A^ch x 8 A'ch unitary scatterings matrices 
Sa, Sh, and Sn, respectively. 



*2 



= Sa 



= 5, 



N 



= Sh 



*3^ 



(19) 



Each of the matrices Si {i — a, b, N) is further decom- 
posed into 4A'chx4A"ch reflection (r^, r'^) and transmission 
{ti, t[) matrices. 



S,^ 






(20) 



Further, the scattering matrices Si are diagonal in 
electron-hole space. 



S^ie) 



S^ie) 
S*{^e) I ' 



(21) 



whereas, in spin space, Sn is proportional to the 2x2 
identity matrix, Sb = d^a%(Sh\ , Sh\) is diagonal, and Sa 
reads 



e -tZvl ( Sa\ \ .Zvl 



(22) 



Finally, the full scattering matrix S of the trilayer can be 
calculated by combining Sa, Sn, and St using Eqs. ( |l9| ) 
and (Efl). One thus finds 






S 






5 



r i' 
t r' 



(23) 



with, 



i = tb(l-r^rt)-iijv (24) 

X [l - r'j'^nil - r'^ny^tpf - r'^r„] ta, 

r = ra +t'^[rN +rb{l- r'ffrby^]tN (25) 

X [l - r'^{rN + nil - r'j^n)-^)] '^ ta, 

and similar expressions for r' and t' . 

The scattering matrices 5^, S^, and Sb are the input 
parameters of our approach. While valuable insight can 
be obtained by using simple ansatzes for these matrices, 
detailed knowledge of the scattering matrices jfor FN in- 
terfaces is available from ab-initio calculationsEj and an- 
alytical results are known for the statistical distributiaa. 
of Sn for the case of a disordered normal metal spacer .tJ 

At the interfaces with the superconducting contacts, 
electrons incoming from the trilayer are reflected as holes, 
and vice versa. This process, known as Andreev reflec- 
tion, is described by the 8Nch dimensional scattering ma- 
trix Sa, 






-qjL 



Sa^ 



rAil) 





rAi-l) 



(26) 



In the limit that the superconducting gap Ao is much 
smaller than the Fermi energy Ep, ta reads, in electron- 
hole grading. 



rA{4>) = a(e) 







icFye 




(27) 



where a(e) = e-*^'-<=™^(<^/'^o) = c/Aq - i^Jl - e^/A^ and 
tjy acts in spin space. The Andreev reflection matrix rA 
is diagonal in the space of transverse modes. 

From Eq. (^) and Eq. ( P6| ) one finds the equation 



SaS 









(28) 



from which it follows that the spectrum of the Josephson 
junction is given by the solutions of 



det[l-5A(e)5(e)] = 0. 



(29) 



Following Ref. E7l one can then express the free energy 
F in terms of Sa and S. Taking derivatives to the su- 
perconducting phase difference and to the angle 9 one 
arrives at expressions for the Josephson current / and the 
equilibrium torque t°'^^ , 



2e ^ , 
/ = —kT > Tr 

n=0 



V^ rr. W ^^^■n)S(iuJn ) 



SA{i^n)S{iu:n)'' 



^cqu ^ ^rj, y- ^^ SA{l^n)%[l^n) 



n=0 



l-5A(it^„)5(«w„)' 



(30) 
(31) 



Here k is the Boltzmann constant, T the temperature, 
and iOn — (2n -I- l)TTkT are the Matsubara frequencies, 
n — 0, 1, 2, . . .. In the next section, Eqs. ( |30| ) and ( |3l| ) 
are the starting points for our calculations of the super- 
current and magnetic exchange interaction. 



IV. JOSEPHSON INDUCED TORQUE 

In this section, we discuss the Josephson current and 
Josephson-effect induced torque for various simple mod- 
els for the ferromagnetic layers Fa and Fb and the nor- 
mal spacer layer N. We consider the case of a short 
Josephson junction, i.e., we suppose that the length of 
the trilayer is smaller than the superconducting coher- 
ence length, or, equivalently, that the inverse dwell time 
inside the trilayer (the "Thouless energy" ) is larger than 
the superconducting gap Aq. For a short junction, we 
can neglect the energy dependence of the scattering ma- 
trices Sa, Sb, and Sn, and evaluate them at the Fermi 
level Ep- 



A. Toy model 

To illustrate the origin of the Josephson effect induced 
magnetic exchange interaction, we first describe a simple 
(toy) model for the scattering properties of Fa, Fb, and 
N: We assume that both majority and minority electrons 
are transmitted perfectly through the two ferromagnetic 
layer Fa (Fb), but pick up phase shifts that differ by an 
amount (3a (Pb) as a result of the Zeeman coupling to the 
exchange field inside the magnetic layer. Transmission 
through the normal metal spacer layer is ballistic as well. 
Such an assumption corresponds to a WKB treatment of 
the exchange field for the case where the exchange field 
does not depend on the transverse direction. The phase 
shifts depend on the transverse mode. In terms of the 
potential V{f) of Eq. (||), they are given by. 



f3a 



(b) 




\/E — Vniaj — \J E — V^nir 



(32) 



where E = Ep — E± is the longitudinal component of the 
kinetic energy (which depends on the mode index). This 
model is simple enough so that one can solve Eq. ( |29| ) 
directly (see, e.g., Ref. g8|). Up to terms independent of 
(j), the Free energy is then given by 



F 



Nci^kT i 



log cosh 



Aocos[((^±7)/2] \ 
2kT / 



(33a) 



where we abbreviated 



7 = arccos f cos(/3a + (3b) cos^ - + cos{(3a - Pb) sin^ , , 

(33b) 

and where (. . .)ch indicated an average over the trans- 
verse modes. Equation (33) reduces to the result of Ref. 
ps] for the case of a ballistic point contact. 

As an illustration, the contribution to the zero tem- 
perature free energy F from one transverse mode is plot- 
ted in Fig. ^ for a generic choice of the phases f3a and 
Pb- We note that F indeeds depends on both (j) and 9. 
While the ^-dependence of F leads to the existence of 
an equilibrium charge current / through the junction — 
the Josephson current — , the 6'-dependencc of F causes 
an equilibrium spin current between Fa and Fb, i.e., a 
magnetic exchange interaction. We also note that the cj) 
and 6 dependencies of F are of comparable size, which 
allows us to estimate the equilibrium exchange interac- 
tion as t'^'^" ~ hlcrit/2e, where /dit is the critical current 
of the Josephson junction. 

The observation that the (j> and 6'-dependences of F are 
of comparable magnitude is valid for arbitrary Pa and Pb- 
However, the location of the minima and maxima in F 
depends on whether 7 is smaller or larger than 7r/2: The 
minimum of F at fixed 9 shifts from (j) = to (f) = tt 
if Pa and Pb (and 9) are such that 7 exceeds 7r/2. A 
minimum of F for (j) = tt corresponds to a tt junction. 
Since 7 interpolates from Pa + Pb to Pa — Pb as one in- 
creases the angle 9 between the two magnetic moments, 
the TT-junction behavior can be induced by rotating one 
magnetic moment with respect to the other if the phases 
Pa and Pb are sufficiently large. Similarly, by varying 
(f>, one can switch the minimum of the free energy from 
9 — to 9 — TT, thus favoring a parallel or antiparallel 
configuration of the magnetic moments. 

Figure represents the contribution from only one 
transverse channel. As different transverse channels have 
different phase shifts Pa and Pb, their contributions to 
the supercurrent and to the magnetic exchange interac- 
tion do not need to add up constructively. In order to ob- 
serve an appreciable supercurrent and/or a supercurrent- 
induced magnetic exchange interaction, the ferromag- 
netic layers must be sufficiently thin or the exchange field 
must be sufficiently weak that the phases Pa and Pb typ- 
ically do not exceed unity — so that all contributions to 



/ or T^"^^ add up constructively. We wish to point out 
that this is not an impossible condition to meet. In fact, 
the same condition applies to the existence of a super- 
current through a magnetic Josephson junction with a 
single ferromagnetic layer. Such supercurrents have been 
observed, see Ref. |lj Moreover, it is important to realize 
that it is only the difference of phases between majority 
and minority electrons that plays a role. Any common 
phases are cancelled out as a result of the Andreev scat- 
tering. As a result, the magnitude and sign of the Joseph- 
son current and the exchange interaction do not depend 
on the phase shifts picked up in the normal-metal spacer 
layer. This is very different from the standard magnetic 
exchange interaction, where the sign and magnitude of 
the interaction depends sensitively on the thickness of 
the normal-metal spacer layer. 




180 



FIG. 3. Free energy _F as a function of 9 and 4> at zero tem- 
perature for the toy model |3^, with fja = 7r/3 and Pb = n/8. 
F is plotted in unit of A'^chAo. 



B. Chaotic normal layer 

Let us now turn to a more realistic model where re- 
flection processes occurring inside the ferromagnetic and 
normal-metal layers and at their interfaces are fully taken 
into account, and where scattering from impurities causes 
the transverse modes to be mixed. In particular, we want 
to study the effect of spin filtering: the fact that majority 
and minority electrons have different transmission prob- 
abilities for transmission through a ferromagnet layer. 
(Spin filtering is the dominant source of nonequilibrium 
current-induced torque for FNF trilayers with normal- 
metal contacts.) The Josephson current is expected to 
decrease with increasing reflection inside the FNF tri- 
layer, and with an increasing amount of spin filtering. 
(Since the Josephson current is carried by Cooper pairs, 
both the minority and the majority electrons must be 
transmitted in order to get a current; see Ref. ^ for a 
discussion of this effect in the context of an FNF trilaycr 
with one superconducting contact.) 

Here we assume that the scattering matrix of the nor- 
mal layer is drawn from the circular orthogonal ensemble 



from random matrix theory, i.e., in particle/hole grading 
we set 



Sn = diag (So ® I2, 5*0 ® Is), 



where I2 is the 2x2 identity matrix in spin space and Sq 
is a 2A^ch X 2iVch unitary symmetric matrix chosen with 
uniform probabihty from the manifold of 2Nch x 2A*'ch 
unitary symmetric matrices. The ensemble of scatter- 
ing matrices corresponds to an ensemble of FNF trilayers 
with the same Fa and Ff,, but different disorder configu- 
rations in N. The circular ensemble is appropriate for a 
trilayer where the normal part would be, fop-example, a 
dirty metal grain or an amorphous material.Ej When the 
transmission probabilities of the ferromagnetic layers are 
small, the circular ensemble can be used for an arbitrary 
diffusive spacer layer .c3 Our choice of a symmetric scat- 
tering matrix implies that the amount of magnetic field 
leaking into the normal layer from the ferromagnets Fa 
and Fb must be sufficiently small, so that time-reversal 
symmetry is preserved inside N. (If time-reversal sym- 
metry is fully broken in iV, both the supercurrent and 
the Josephson-effect induced torque will vanish to lead- 
ing order in Nch-) 

We are interested in the limit of a large number of 
channels iVch- (For metals, A^ch is already of the order of 
lO'^ for contacts with a width of a few nm.) For large iVch, 
sample-to-sample fluctuations of the current and torque 
are much smaller than the ensemble average, so that the 
ensemble averaged current / or torque -f°i" is sufficient 
to characterize a single sample. In order to calculate the 
ensemble average, we first rewrite Eq. ( P9|) in a slightly 
different form. 



det[l-5fs(e)^iv]=0. 



(35) 



where Spsi^) is the 8iVch x SA'ch matrix describing the 
combined effect of scattering from both ferromagnetic 
layers backed by the superconductors, as seen from the 
normal spacer layer. The matrix Sps has spin, parti- 
cle/hole {e — h), channel, and ''a — b" degrees of freedom, 
the latter grading referring to whether scattering is from 
Fa or from Ff,. In the a — b grading, Spsi^) reads. 



Sfs = 



f Sps.a 
^ Sfs, 



with 



Sfs. a =ra+ta — 



rA{(l)/2)ra 



(36) 



rA{H2)C (37) 



and 



Sfsm ^ n + t'ij 



1 - rA{-(f>/2)n 



rA{-(jy/2)h. (38) 



With these notations, the ensemble averaged Josephson 
current reads, 



n=0 ^ 



(34) with 



Q = 



1 



1 — SfsS_ 



-Sfs- 



N 



(39) 



(40) 



The average Q is computed in the appendix using the 
method of Ref. ^. The results of that calculation is a 
self-consistent equation for t/, analogous to the Dyson 
equation for the average Green function in a standard 
impurity average. 



1 



1 — 5fsS 



S 



FS 



where 



S = 



1 



1 + 4[P(5)]2-1 



(41a) 



(41b) 



2V{g) 

and P is a projection operator. In e — /i space, P reads 



Ahe Ahh 



— tr^.„a6 ( ? 



2^ch 

® lab ® ljV,h 



Aeh 





(41c) 



where the trace tT: N^^^ab ... is taken in channel and a — b 
space, but not in spin space or particle/hole space. 

Equation ( p4| ) reduces to a self-consistent equation for 
the 4x4 matrix E. This equation remains fairly com- 
plicated and, in general, has to be solved numerically, 
even when Sfs is diagonal in channel space. In the limit 
where Sfs is close to the identity matrix, i.e., when both 
ferromagnetic layers are poorly transparent and reflect 
majority and minority electrons with almost the same 
reflection phase, a further simplification of Eq. ( [4l| ) is 
possible. Assuming that Sfs is diagonal in channel space 
(i.e., the ferromagnetic layers do not mix channels), ex- 
panding Sfs = 1 + 5Sfs + 0{5SFsY^ ^"^^ defining the 
4x4 matrix X = l/{2Nch)tT:N^i,.ab5SFSi Eq- (0) reduces 
to 

Xhe + ^heXee + Xhh^he + ^heXeh^he — 0, 

Xeh + 'SehXhh + Xee^eh + 'SehXhe'Seh = 0, (42) 

while Eee = ^hh = since E = P(E). Equation (|4|) 
shows that for opaque FN interfaces and for diffusive 
scattering from the normal spacer layer, there is only 
a restricted number of parameters (i.e., the free parame- 
ters of Xeh ; Xhe is related to Xeh by particle- hole sym- 
metry) that determines the supercurrent and the equilib- 
rium torque. 

We could only obtain a solution in closed form in the 
case where transmission and reflection amplitudes of the 
ferromagnetic layers were real, i.e., still allowing differ- 
ent transmission and reflection probabilities for majority 
and minority electrons (spin filtering), but without spin- 
dependent phase shifts in Fa and Ff,. Introducing the 



mode-averaged transmission probabilities Tai {Tai) for 
majority (minority) electrons, 



the spin-averaged transmission probability Ga ~ {Ta^ + 
Tai)/2, and the geometric mean 7a = {Ta^TaiY^'^, with 
similar definitions for Tbi, Tf,|, Gt, and 7f,, we find 






VA^+^« 
7q76 



(43) 



V7a + 7h + 27a76 COS (/) -K (w„/Ao)2(Ga + Gb)^ ■ 

Equation ( [43|) is the generalization of Equation (24) of 
Ref. |2^ to the case of contacts with spin-dependent trans- 
mission, in the short junction limit (Thouless energy 
much larger than superconducting gap Aq). In the limit 
of zero temperature, it reduces to a complete elliptic inte- 
gral of the first kind. We note that the Josephson current 
( [43| ) does not depend on the angle 9, so that there is no 
equilibrium torque in this case. Since the assumption 
underlying Eq. (|4^ ) was that the electrons do not pick 
up phase shifts in the ferromagnetic regions, we conclude 
that spin- filtering alone (i.e., the fact that majority and 
minority electrons have different transmission probabil- 
ities) is not enough to create a (Josephson current in- 
duced) equilibrium magnetic exchange interaction. We 
numerically checked that this conclusion still holds when 
Sps is not close to unity. 



C. Numerical results 

For a numerical solution of Eq. ^ that accounts for 
the fact that majority and minority electrons experience 
different phase shifts while scattering from the ferromag- 
netic layers, it would be desirable to have detailed knowl- 
edge of the scattering matrices Sa and Sb of the ferromag- 
netic layers. These scattering matrices can, in principle, 
be calculated from ab-initio calculations, see, e.g., Ref. 
Pq . However, complete data for all amplitudes (phase 
shifts and probabilities) are not available in the litera- 
ture. Therefore, we choose an ansatz for Sa and Sb that 
is close in spirit to the toy model of Sec. IV A. The use 
of a simple ansatz can be partly justified by Eq. (42), 
which shows that it is only a finite number of parameters 
that determines the supercurrent and exchange interac- 
tion, not the entire matrix Sa or Sb- For our ansatz, 
we assume that these scattering matrices are diagonal in 
channel space and we neglect any channel dependence. 
We further assume that most of the reflection processes 
take place at the FN interface. Without loss of general- 
ity, we may set the phase picked up by minority electrons 
while traversing the ferromagnetic layers equal to zero. 
Then the difference between minority and majority elec- 
trons is fully described by the phases f3a{b) picked up by 



majority electrons. This leads us to the ansatz (in left- 
mover/right-mover space) 






LAfci 



s, 



fci.ee 



_ ( \J^ - Tbi iy/Tbj 



v/i^Ti 



6i 



•^IWch: 



(44) 



and similar equations for S'at.ee and Sai.ee, while the scat- 
tering matrices Sb^.hh, Sbi,hh, 'Sa.yih, Saijih are given by 
the complex conjugates, cf. Eq. (pi|). 

With this model for the scattering matrices Sa and Sb, 
a typical plot of the supercurrent at (j) — tt/2 as a function 
of 6 is shown in Fig. El Wc have taken the values of the 
parameters Tqi, Tai, Tb^; lifvi Tbi from realistic estimates 
for a Co-Cu-Co trilayer,t30 while we fixed the phases Pa 
and (3b arbitrarily. Although the choice of parameters is 
specific, the observation that the Josephson effect induces 
a magnetic exchange interaction between the ferromag- 
netic layers was found to hold for any generic choice of 
scattering parameters. [The only exception being the 
case discussed around Eq. (43), for which all scattering 
phase shifts are either or it.] Further, we found that 
when the phase difference between minority and major- 
ity electrons becomes of order unity, the variation of the 
supercurrent / with 6 is of the order of the critical cur- 
rent, so that, up to a numerical factor, the magnitudes 
of maximum equilibrium torque and critical current are 
related as r'"'" ~ ?i/crit/2e. 




-180 



180 



FIG. 4. Supercurrent / at phase difference (p = 7r/2 be- 
tween the superconducting orderparameters, as a function of 
the angle between the moments of Fa and Fb- The cur- 
rents is measured in units of 2eA'chAo/7i. We set Ta]- = 0.68, 
Tai = 0.29, Tb-f = 0.68, Tti = 0.29, following ab-initio stud- 
ies of Ref. Eq for a Co-Cu-Co trilayer. The phases f3a and 
/3b were arbitrary set at f3a ~ 357r/180 and f3b = 1157r/180. 
A:r = 0.1Ao. 

We now turn to a slightly different model for the nor- 
mal layer which we study by doing the disosder average 
numerically using Eq. (BQ). In this model,E3 which was 
also used in Refs. |^ and|l% Sn is given by Eq. (34), where 
the 2A'ch X 2iVch scattering matrix Sq is parameterized, 
in a-b grading, as 



So 



u 

f/T 



(45) 



+ J: 



sin 261), (46) 



Here U is an N^h x Net unitary matrix, uniformly dis- 
tributed in the group of unitary Nch x Nch matrices. This 
model, in which the normal metal spacer mixes the trans- 
verse modes, but does not cause any backscattering is 
appropriate, e.g., for rough FN interfaces. As the same 
model was considered for quantitative estimates in Refs. 
@ and E2I this choice for Sn can be used for a quanti- 
tative comparison of the Josephson-effect induced equil- 
brium torque and the nonequilibrium torques considered 
in Refs. g and |l2|. Results are shown in Fig. ^ for the 
same choice of parameters as in Fig. ^ 

For all choices of Sa and Sb that we considered, we find 
that the results are well described by the phenomenolog- 
ical relation 

9t„^i" h dl 

itr = ^^- ^ch Ao sin Ji sin 9 

where the constants Ji and J2, which are analogous to 
the quadratic and bi-quadratic coupling constants in the 
standard magnetic exchange interaction, depend on the 
detailed choices for Sa and Sb- Several properties of this 
phenomenological relation are worth while mentioning, 
(i) The torque induced by the Josephson current is pro- 
portional to the number of transverse channels iVch, and 
hence to the width of the trilaycr A^ch- (This property 
holds if the ferromagnets are sufhciently weak or thin, so 
that the phase difference experienced by majority and mi- 
nority spins is ^ 1, see the discussion at the end of Sec. 
IV A.) This should be contrasted with the regular ex- 
change interaction which does not increase with increas- 
ing A"ch for a disordered normal-metal spacer .El Thus, for 
wide junctions, the Josephson torque is parameterically 
larger than the standard magnetic exchange interaction, 
(ii) The torque r*"^" vanishes &t 9 — Q and 9 — ir, ir- 
respective of (j), since for these angles only spin currents 
parallel to the magnetic moments play a role, which are 
conserved, (iii) Similarly, the Josephson current / van- 
ishes at ^ = and (p = tt, irrespective of 9. (iv) The 
Josephson current depends on the angle 9. In the case 
shown in Fig. ra, the junction shows 7r-junction behavior 
for 9 — 0, which disappears when 9 approaches tt. The 
relative strength of the 6'-dependent and ^-independent 
part of the current varies with the phases Pa and Pb, and 
the TT-junction behavior is not necessarily there for all 
choices of Pa and Pb- 



V. CONCLUSION AND DISCUSSION 

We have shown that for a ferromagnet-normal-metal- 
ferromagnet (FNF) trilayer coupled to two superconduct- 
ing contacts, the Josephson effect enhances and controls 
the magnetic exchange interaction between the magnetic 
moments of the two ferromagnetic layers. 
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FIG. 5. Josephson current / and equilibrium torque 
T°''" for the model mS) for the scattering matrix of the nor- 
mal-metal spacer layer, and with the same choice of Sa and 
Sb as in Fig. kl From up to down, the panels contain: 7(0) for 



= 7r/2, I{9) for (j) = 7r/2, r|'i"((?!>) for 6 = 7r/2 and 



\e) 



for (j) — -K 12. For the choice of parameters used for this figure, 
the Josephson current and the equilibrium torque can be fit- 
ted with the phenomenological relation (|46[), with J\ = 0.007 
and J2 = 0.025. Torques are given in units of AoA'^ch and 
currents in units of 2eAoA'ch/?i. ^T = O.IAq. 



The Josephson-effect induced torque bears important 
similarities and differences with the nonequilibrium and 
equilibrium torques in an FNF trilayer with two normal 
metal contacts: 

• The Josephson-effect induced torque points perpen- 
dicular to the plane spanned by the directions rria 
and rhb of the magnetic moments of the ferromag- 
netic layers, like the standard equilibrium exchange 
interaction. On the other hand, the nonequilibrium 
torque mainly lies inside the plane spanned by fha 
and rhb- 

• For the Josephson-effect induced torque (or the 
standard equilibrium exchange interaction) to ex- 
ist, transmission through the FNF junction needs 
to be phase coherent. Moreover, existence of the 
Josephson-effect induced torque requires that ma- 
jority and minority electrons experience different 
phase shifts upon transmission through or reflec- 
tion from the ferromagnetic layers. The nonequi- 
librium torque, in contrast, only needs spin filtering 
(different transmission or reflection probabilities for 
majority and mipxidty electrons), while coherence 
is not important .H'Lj 

• Like the nonequilibrium torque, the Josephson- 
effect induced torque is carried by states close to 
the Fermi energy. The standard equilibrium torque 



has contributions from states throughout the con- 
duction band. 

The equihbrium torques f^'^" and f^'^" on the mo- 
ments of both ferromagnetic layers Fa and Ff,, re- 
spectivly, are equal in magnitude, but opposite in 
direction, f^'i" 
the nonequilibriuni torque 



f^'^". No such relation holds for 



• The sign and size of the nonequilibrium torque 
is controlled by the direction of the current. In 
contrast, the sign of the Josephson effect induced 
torque is set by the superconducting phase differ- 
ence (j) and by the details of the scattering phase 
shifts from the ferromagnetic layers; it is not re- 
lated to the direction of the supercurrent in any 
direct way. However, the order of magnitude of the 
Josephson effect induced torque is set by the size 
of the critical supercurrent, r'^'^" ^^ hlcrit/2e. 

We close with a discussion on the effect of the Joseph- 
son induced exchange interaction on the dynamics of 
the magnetic moments. Typically, one of the two fer- 
romagnetic moments (say the moment rfii, of the layer 
Fb) is fixed by anisotropy forces, and the torque is stud- 
ied through its effect on the moment rha of the "free" 
layer Fa . The usual method to describe the dynamics rfia 
in the presence of the current induced torque is viapthe. 
phenomenological Landau-Lifshitz-Gilbert equationE'tj 
The result of such a calculation is a critical value of the 
torque necessary for switching rua with respect to the 
fixed moment m^. This program was carried out for the 
non-equilibrium torque acting on a trilayer connected to 
two normal electrodes in Refs. ^ and ^. The effect of 
the equilibrium torque considered here is simpler, as it 
admits a formulation in terms of the total energy of the 
system. Using the phenomenological relation (Eq), the 
Josephson induced exchange interaction corresponds to 
an energy gain Sf per unit area equal to 



df = -^C08(j) [Jirha ■ rhb + J2{rha ■ rhbf 



(47) 



The criteria for switching the orientation of rfia is that 
the magnitude of this interaction energy exceeds that of 
the work done against anisotropy forces acting on each 
layer (arising from shape, crystalline structure, etc.). For 
a 2 nm thick Cobalt layer, those are of the order of 10""^ 
Jm^^ and can be decreased by up to two orders of mag- 
nitude for if Cobalt is replaced by PermalloyEJ NigiFcig. 
On the other hand, for Aq ~ lOK (as is the case for 
Niobium), Xp ~ lA and Ji « 0.01 (the value found in 
our toy model simulations with slightly optimized values 
for the phase shift differences Pa and Pb), the Josephson 
induced interaction is of the order of 10~* Jm~^. Hence, 
we estimate that control of the relative orientation of rfia 
and rhb should be experimentally accessible provided the 
local anisotropy forces are kept at a minimum. 

When the anisotropy forces are so small that switch- 
ing of the magnetic moments becomes a possibility, the 



ac Josephson effect should provide a clear signature of 
the switching of the ferromagnetic moments through the 
sensitivity of the supercurrent to the angle 9 between 
rha and rhb- A possible scenario is sketched in Fig. ]^: 
the equilibrium current observed in the Josephson effect 
should exhibit periodic shifts when the switchings oc- 
cur. In order to evaluate the fastest time scale at which 
the switching of the moments can occur, we return to 
the Landau-Lifshitz-Gilbert equation for the dynamics 
of the magnetic moments. As before, we suppose that 
rhb is kept fixed by a strong local anisotropy field, while 
the anisotropy field acting on rha is negligible. Neglect- 
ing the bi-quadratic oapBling J2 , the dynamics of rha can 
then be described by,BE3 



drha 



rUa X 



jHjirha ■ rhb)rhb - a 



drha 



dt 



(48) 



Here 7 = gfJ-s/f^ is the gyromagnetic ratio, a the Gilbert 
damping coefhcient and Hj is the effective exchange field 
representing the Josephson-effect induced torque, 



Hj 



•h 



ma 



■COS0(t), 



(49) 



where da is the thickness of the layer Fa- Neglecting the 
time dependance of the superconducting phase difference 
(/), this equation is easily solved, 



tan0 = e 



-Tt 



tan^n 



(50) 



with r = ^Hja/{\ + q;^). The switching rate T is max- 
imum for a — 1 where, using the numerical values con- 
sidered above and \ma\ = 1.6 lO^A/m (for cobalt), we 
find r w 10 GHz. On the other hand, typical voltages 
used in ac Josephson experiments are of the order of a 
few /iV, which corresponds to frequencies of the order 
of a GHz. Therefore, for these frequencies, it should be 
possible, in principle, to observe the periodic switching 
of the magnetization orientation as suggested in Fig. O. 




eVt 



FIG. 6. Sketch of the current as a function of time when 
a small voltage V is applied across the junctions. When 
the Josephson effect induced interaction exceeds the local 
anisotropy field, the relative configurations of the two mo- 
ments switches. 
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APPENDIX: DERIVATION OF EQ. @ 



In this appendix, we calculate the ensemble average Q 
where Q is defined in Eq. (^). The ensemble average 
amounts to the calculation of an average over the circu- 
lar orthogonal ensemble from random matrix theory (the 
manifold of unitary symmetric matrices), 



_ (a) 



/ 



G — dSo Sfs 



1 



1 — SnSfs ' 



(Al) 



where, in particle/hole (e — h) grading, Sn = diag {Sq (^ 
I2, 'S'o (81 I2), I2 being the 2x2 identity matrix in spin 
space and S'o being an 2Nch x 2A^ch symmetric unitary 
matrix, and dSo is the invariant measure for integration 
over the circular orthogonal ensemble. To perform the 
average over So, we use the diagrammatic technique of 
Ref. ^ and calculate G to leading order in 1/iVch- First, 
G is expanded in powers of Sps , 



G — Sfs + SfsSnSfs + SfsSnSfsSnSfs 



(A2) 



and the corresponding terms are associated with dia- 
grams: a full line corresponds to a factor Sfs and a dot- 
ted line to factor Sn, see Fig. 0a. At each dot one sums 
over a latin index ranging from 1 to 2Nch, representing 
the channel space, and a greek index ranging from 1 to 
4, representing the combination oi e — h and spin space. 
According to the diagrammatic rules of Ref. ^, the av- 
erage is then done by connecting all dots by thin lines, 
summing over all possible ways of pairing up the dots. To 
leading order in 1/Nch, only planar diagrams contribute, 
i.e., the diagrams where the thin lines do not cross. When 
two dots are connected, the corresponding latin indices 
are identified and summed over, and the constraint is im- 
posed that two greek indices involved have to be different 
in e—/i space (i.e., if one index corresponds to e, the other 
one has to correspond to h). The rationale for this con- 
straint is that only contractions that involve S'o and its 
complex conjugate are allowed in the average. Finally, a 
weight factor is associated with each diagram: each cycle 
formed by an alternation of dotted lines and thin solid 
lines in the diagram contributes a factor Wi, where i is 
equal to half the number of dotted lines contained in the 
cycle. The Wi are tabulated in Ref. p9| for the purpose 
of this integral, we only need their generating functionE^I 






w,y 



-1 = (v/(27V,h)2+4z - (2iVeh)) /2z. (A3) 



The weight factors Wi with i > 1 correspond to contribu- 
tions to the average that go beyond a Gaussian evaluation 
using Wick's theorem. 

The first non vanishing diagrams are shown in Fig. [7|b. 
They correspond to 



G = Sfs + 2NchWi Sfs P{Sfs) Sfs 



(A4) 



G = — + —•— + _ ..- ^. .— + 

ij iklj iklmnj 
apa^^pa)J.|J.ppi3 



G = 



+ W. 



+ 



(b) 



ij ikmmkj 




(c) 



G 



(d) 




P(G) 



P(G) 



P(G) 



FIG. 7. (a) Diagrammatic representation of Eq. (A2). 
(b) The first two diagrams contribution to Q. (c) General 
structure of diagrams contributing to Q. (d) General structure 
of diagrams contributing to E. 



The projector operator P was introduced in Eq. ( 41c ). 
It implements the constraint that only greek indices rep- 
resenting e and h degrees of freedom can be contracted. 
Now, we are ready to obtain the self-consistent equa- 
tion ( [4l| ) for G. Except from the zeroth order (first term 
in Eq. (A4), all the diagrams involved in G have the 
structure shown in Fig. |7|c, where the boxes stand for all 
possible allowed contractions. The important point here 
is that there is no thin line connecting the two differ- 
ent boxes since we are only considering planar diagrams. 
The sum over all the different contractions represented by 
the right box gives all the possible diagrams and there- 
fore equals G itself. Denoting the left part of Fig. ^ by 
E, we have 



G = Sfs + Sfs^G, 



(A5) 



which is Eq. (41a). The diagrams contributing to S have 
the structure shown in Fig. Md. They contain 2n— 1 build- 
ing blocks, n ~ 1,2, . . ., with weight factor Wn- (Only 
odd numbers appear, because of the constraint that only 
indices belonging to So and its complex conjugate can be 
contracted.) Each building block can be identified with 
P(^). Hence, 



S = ^I^„Tr [2N,hPiG)Y 



(A6) 



n=l 
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which leads to Eq. (41b) if we use the generating function 
( [A^ ) for the weight factors Wn- 
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